Numerical solutions to the cosmological 3-fluid problem 
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Abstract 

We show that, for the scalar field cosmology with exponential potential, the set of values of the 
coupling parameter for which the solutions undergo a transient period of acceleration is much larger 
than the set discussed in the literature. The gradual inclusion of ordinary and dark matters results in an 
everywhere, but near the origin, smoother and right shifted (along the time axis) acceleration curve. 
For the 3-fluid problem, the energy density need not exhibit a plateau during the acceleration period. 
Much excess in the dark matter and/or ordinary matter energy densities would lead the universe to 
undergo an eternal deceleration expansion. 

PACS numbers: 98.80.-k, 95.35.+d, 95.36.+X, 04.25.D- 



The problem. To account for the present period of acceleration of the universe many models have been 
^ ■ suggested so far HI- ifTSl . In this work we are concerned with the models corresponding to spatially 

flat Friedmann-Roberston-Walker (FRW) universes where the source is (1) a homogeneous scalar field 
^ ' coupled to exponential potential [representing dark energy (DE)] |l8][l4l and references therein or (2) 



a homogeneous scalar field coupled to exponential potential along with the dark matter (DM), which is 
modeled by a dust fluid, and a perfect fluid representing ordinary matter (OM) lITSl and references therein. 
^ I This is the so-called 3-fluid problem. 

' The problem of a scalar field coupled to exponential potential has been solved in (i-dimensional by 

decoupling the system of equations governing the dynamics [8]. If the potential is of the form V = Vqc^^'^ 
where (p is the scalar field, the solutions have been classified according to the coupling constant A fl4l . 
^ ■ The problem with a piecewise exponential potential in 4-dimensional has also been solved 1 14|. 

^ . The 3-fluid problem that involves a scalar field (DE), a dust fluid (DM) and a perfect fluid has been 

treated analytically in 1 15 1. The method developed in |[T5l and references therein that is based on sym- 
metry assumptions has led, as was shown later in fr6l, to a flawed solution to the 3-fluid problem, in that 
the solution derived in [ 15 1 does not satisfy the field equations. 

Thus, to our knowledge both problems of deriving analytic solutions or developing new methods 
to tackle to 3-fluid problem remain open. As we shall see below, even the linear case is not easy to 
handle. The purpose of this work is to provide numerical solutions to the 3-fluid problem with exponential 
potential and compare them to the solutions of the scalar field problem with similar potential. 

To write explicitly the solutions for the normal (not phantom) scalar field problem we use a slightly 
different notation but the same conventions as in |[T4l with the exception of taking V = Vbe^'*''^ (instead 
of y = Voe+^'t'). For the case of the scalar field alone we use the subscript s to denote any variable. For 
instance, we write = Vbe^'*''^' and = for the scalar field problem and V = Voe^^*^ and = p for 
the 3-fluid problem. 

Consider the flat (k = 0) FRW cosmological model for the universe equa- 
tions governing the motion of {as,(j)s) are the FRW equation HJ = Ps and the Klein-Gordon one, where 



1 



Ps is the energy density and Hg = ds/us is the Hubble variable (i = dx/dt) 



• 2 A 2 

a? ~ 2 



^ = ^+Voe-^«'^ (1) 



+ =0. (2) 

a* 

The solutions are called hyperbolic if A < 3\/2 or trigonometric if A > 3\/2 ||T45. If we let as = e^'/^, 
0, = V2Ws/3, a = Vl8-A2/(3\/2), j8 = 7^2-18/(3^2), and £± = 6/(6± \/2A) the hyperbolic and 
trigonometric solutions are given by, respectively 

Us = 2T/a + ln(l + me^^ar) ^ jj^^ ^ _ ^e-2aT>) 

Ws = V2Xr/{3a) + ^_ ln(l + me"^"^) - £+ ln( 1 - me^^"^) (4) 

Us = i- ln[cos(j8T)] +41n[sin(/3T)] (5) 

Ws = £_ln[cos(j8T)] -^+ln[sin(j8T)], (6) 

where we have omitted additive constants. Here m is a real constant and T is the new time coordinate 
defined by 

t = ^^e-^^«''/^ (7) 

If m 7^ 0, then m may be set equal to ±1 by shifting the values of T. For all A < 3\/2, it is straightforward 
to show that the range of T extends from a finite value (usually taken 0) to infinity, and for A < 3\/2, 
0< T< 71/(2/3). 

The 3-fluid problem is also governed by the Klein-Gordon equation ^ (with the subscript s removed) 
and the FRW equation = p with H = d/a and p = p^+ Pdm + Pom is the sum of the energy densities 
of the three fluids (scalar field, dark and ordinary matters). If the scalar field is normal, its energy density 
is given by 

P0 = y+Voe-^^. (8) 
The DM is treated as a dust (without pressure) while the OM is a perfect fluid with an equation of state 

PoM = {y-l)poM (1<7<2). (9) 

Using the three conservation equations (for each fluid T^^ -y = where T^^ is the stress-energy tensor) 
and Q, one arrives at Iil5j 

PoM = Ca"^^, pDM = £'/a^ (10) 

where C > and £ > are integration constants. The 3-fluid problem reduces to that of the scalar field 
if C = £ = 0. The equations of motion take the fonr0 

aa2_ 1^3^2_^^^3g-A0_j^^-3(r-l) (11) 

^ + 3-(^-AVoe-^^ =0. (12) 
a 

Now we set a = eP/^ and = \/2W /3 reducing the system to 

9Ee-^ = U^-W^- Woe-^^^/^ - QCe"^^ (13) 
2W + 2UW - 3V2XVq&-'^^'^ 1^ = 0. (14) 



'The second teim in Eq. (20) of Ref. Il5l should read 
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As noticed earlier, ifC = E = the system (ITTI ) and ([T2l ) reduces to that of the scalar field problem, which 
is decoupled upon introducing the new variables (m^ , v.,) such that Us = v., + Ug and Wg = Vs — Us and using 
the time T defined in ([7]) [8]. It does not seem easy to decouple the system for C 7^ and E y^O, the case 
to which we are interested. 

For the purpose of the numerical analysis we are aiming to perform, it is not convenient to introduce 
a new time coordinate T defined by t = ^^/V(^e^^^^^^. Rather, we will use the same time T as defined 
in This will allow us to compare solutions to both problems with the same initial conditions. The 
system ([T3l ) and ([T4l) takes the form (with y = dx/dr) 

4£g-f/+v/2Aw,/3 = VoU'^ - VoW'^ - 4Voe-^^(^-"'')/3 - 4Ce-^^+^^^'/3 (15) 

6W" - VlXW^W + 6U'W' - 4v^Ae-^^('^-'^')/3. (16) 

It is possible by a perturbation approach to decouple the system (ITST ) and ([T6l ) if we assume that 
the independent constants C and E are small enough to allow for such an approach. In fact, writing 

U = Us + UcC + UeE -\ andW = Ws + WcC + WeE H it is possible to determine the four time 

functions iUc,UE,Wc,WE), however, their expressions are sizable and a numerical analysis is needed 
anyway. We would prefer applying a numerical approach to the full nonlinear system ([TSl l and ([T6l ) where 
we can consider large values of C and E. 

Numerical solutions. For the scalar field problem, it was noticed in Ref. ||8] that the energy density 
exhibits a plateau during the transient period of acceleration (TR4). We verify that this statement remains 
valid for at leasH \/6 < A < 3\/2 and A > 3\/2. However, the existence of a plateau for the energy 
density does not ensure the occurrence of a TPA: We have noticed the existence of solutions, as are the 
cases for A = Vl and X =2/ -v/3, with no TFA but (twice) the kinetic energy vanishes at the point where 
the energy density has a plateau and the potential energy is maximum. For instance, for A = ^/2 a plot 
of the energy density (with a plateau at T = 0.935), potential and kinetic energies is identical to Fig. 1 of 
Ref. [8]. Thus, the clarification provided in [5| as to the origin of the TPA lacks support within the scalar 
field cosmology itself. As we shall see below, in the extended scalar field cosmology, that is the 3-fluid 
cosmology, the energy density need not exhibit a plateau during the TPA as Fig. [T] reveals. 

For all solutions discussed below we choose Vq = 1 and 7=1.5 and we plot in red (analytic) solutions 
to the scalar tield problem and in blue (numerical) solutions to the 3-fluid problem. Our initial conditions 
are [/(tq) = t/.(To), W{to) = W.(to) and IV'(to) = w;(to) with Tq = IQ-^". 




Figure 1: Case A = 2V2 with m= 1. (a) (left): The energy densities p = /f^ for C = £ = 1 and pj = H^. has a plateau 
around t = 0.65. p exhibits no plateau where the acceleration occurs, (b) (medium): The accelerations a for C = £ = 1 and a.,. 
ds is maximum at T = 0.72. (c) (right): ci for C = E = 100. 

We start by discussing the solution plotted in Fig. 1 of Ref. llH for which A = 2^/2. For C = E = I 
Fig. [11(b) shows plots of dg and d which both exhibit three phases: deceleration-acceleration-deceleration 

^Because of different conventions, the values of (A«, Vor) used in (8) are such that A — y^A^ and Vqs = 6V0. 
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Figure 2: Case 1 = 3 with m = 1. (a) (left): a for C = £ = 1 and a^. (b) (medium): d for C = E = 10. (c) (right): d for 
C=l,E= 100. 



a a , a a a 




Figure 3: Case A = 5. (a) (left): a for C = £ = 1/100 and flj. (b) (medium): ci for C = £ = 1/10. (c) (right): a for C = £ = 1. 

(DAD). For larger values of the constants C = E = 100 as in Fig. [T](c), the TPA is right-shifted and the 
acceleration is smoothed. For much larger C and E the acceleration is completely smoothed, except near 
the origin, and the solution no longer exhibits the DAD phases; rather, it exhibit an eternal deceleration 
expansion [an example is shown in Fig.[3](c) for A = 5.] In Fig. [2] we have similar plots for A = 3. 

These same conclusions extend to the cases A < 3\/2 where no TPA occurs in both problems. In this 
case the acceleration curve, which exhibits deceleration-eternal acceleration (DA) expansions, is right 
and downward shifted with increasing C and E until an eternal deceleration occurs. 

Curiously enough, the trigonometric solutions (jS) and ^ were not discussed deeply in the litera- 
ture |l8][l4l. Contrary to a statement made in llT4l which claims that these solutions exhibit only the DA 
expansions, we have checked that for A > 3\/2 up to (and above) A = 100 these solutions do indeed 
exhibit the DAD expansions, as the red plot in Fig. |3] (a) depicts. Moreover, the acceleration period co- 
incides with a plateau for the energy density. Numerical solutions to the 3-fluid problem are plotted in 
Fig. [3] showing the acceleration for increasing values of the constants. The acceleration phase disappears 
as the constants increase and the solution ends with an eternal deceleration expansion. 

We have checked that the conclusions drawn for the case 7=1.5 extend to the case 7 = 1 .9 and most 
probably to all 1 < 7 < 2. 

Discussion. Solutions to the scalar field problem already predict the DAD expansions for the cosmolog- 
ical evolution of the universe for a wide range of A. Solutions to the 3-fluid problem add new dimensions: 
The parameters C and E along with A allow for fine tuning and horizontal and vertical shifting of the ac- 
celeration and density curves to fit observational data. Excess in the OM density (C large) or/and the 
DM density (E large) results in a smooth acceleration curve, except near the origin, exhibiting no DAD 
phases; rather, the universe would undergo an eternal deceleration expansion. 

No solution with two or more TPA's were found. So, the model, with OM and DM included, still fails 
to predict the early universe inflation. 
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